Re-examining the Verwey transition in Fe3C>4 
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Motivated by recent structural data questioning the adequacy of the charge order (CO)/disorder 
picture for the Verwey transition (at T = TV) in magnetite, we re-investigate this issue within a new 
theoretical picture. Using the state-of-the-art LDA+DMFT method, we show that the non-trivial 
interplay between the B-site octahedral distortions and strong, multi-orbital electronic correlations 
in the half-metallic state is a necessary ingredient for a proper quantitative understanding of the 
physical responses across TV- While weak CO is found to have very small effects on the low- 
T spectral function, the low-T charge gap and the resistivity jump across Ty are quantitatively 
reproduced only upon inclusion of CO in LSDA+DMFT scheme. Our results strongly suggest that 
the Verwey transition is dominantly driven by multi-orbital electronic correlations with associated 
JT distortions on the B-sublattice, and constitutes a non-trivial advance in attempts to understand 
the physics of Fe304. 
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PACS numbers: 71.28+d,71.30+h,72.10-d 

Magnetite is a ferrimagnetic spinel oxide with a high 
Curie temperature, T c = 858 K. Hence, it is viewed as 
an ideal candidate for room temperature spintronic ap- 
plications. Being the oldest magnetic material known to 
humanity, it has nonetheless defied a consistent physical 
understanding to date. It is believed that the Verwey 
transition at Ty = 123 K 0] is associated with an order- 
disorder transition from a charge ordered (CO) state of 
Fe ions on the B sublattice (with Fe 25+5 Fe 25 ~ 5 f or 
Fes ions) to a disordered state (Fe 2 5 at all Fes sites) 
at higher T. One would expect insulating behavior below 
Ty, as indeed observed. However, a consistent picture of 
the electronic state and transport has proved to be very 
controversial. Photoemission (PES) studies lead to con- 
flicting results: Chainani et al. found a finite density of 
states (DOS) at the Fermi level, in agreement with the 
predicted metallic behavior in the high-T cubic phase. 
This is apparently corroborated by the optical conductiv- 
ity studies Q, which reveal clear gap opening below TV, 
and significant spectral weight transfer at higher energies 
at higher TV The "Drude" contribution is estimated to be 
small, with most of the spectral weight spread out over 
the incoherent part. However, Park et al. Q use PES 
and IPES results to conclude that there is no spectral 
weight at Ep above T v . The gap does not collapse across 
Ty, but only shrinks by about 50 meV . This would be 
consistent with a conductivity jump by a factor of 100 
across Ty only if semiconducting behavior is assumed 
above Ty as well. Qualitatively similar conclusions have 
been derived very recently by Schrupp et al. 5], who, 
moreover, use soft X-ray photoemission to determine the 
bulk electronic structure (PES lineshape). Thus, the 
Verwey transition is a semiconductor to semiconductor 
transition, rather than an insulator-metal transition, as 
thought earlier. Other observations like the small carrier 
mobility, the unusual T-dependence of the conductivity, 
and the mid-infrared peak in optics Q suggest electronic 



localisation, but now in both phases. 

The existence and nature of the CO state in Fea04 
is equally controversial. New high resolution powder 
diffraction data find very small charge disproportion- 
ation, —0.1 < 6 < 0.1, while resonant X-ray scattering 
studies even concluded the absence of CO below Ty 0. 
Thus, contrary to traditional expectation, the electro- 
static minimisation required for the CO cannot be the 
dominant mechanism driving the Verwey transition. On 
the other hand, the possibility of orbital order (OO), per- 
haps of a more complex type, has not hitherto been in- 
vestigated. Recall that Fe 2+ is a Jahn- Teller active ion, 
and, in fact, that the low-T distortions of the O octa- 
hedra around B site ions can be viewed as JT distor- 
tions Indeed, a recent LDA+U calculation finds that 
the local occupation of t2 g orbitals on B sites follows the 
octahedral distortion expected from a JT effect. They 
also find a charge disproportionation close to the experi- 
mentally suggested (small) value, but this could probably 
be washed out by dynamical frustration effects beyond 
LDA+U in the real system. Observation of the distorted 
octahedra around B sites in the refined low-T structure 
suggests that orbital correlations and associated JT ef- 
fects in a fully magnetically polarised situation might be 
essential ingredients for understanding the Verwey tran- 
sition. 

An important feature is clear from PES, IPES and op- 
tical measurements across Ty . Characteristic strong cor- 
relation signatures are clearly visible. In particular, the 
detailed forms of the lineshapes in each phase, as well 
as the appreciable transfer of dynamical spectral weight 
over large energy scales across Ty constitute clinching 
evidence for a phase transition between two strongly cor- 
related phases accompanied by strong spectral weight 
transfer (SWT). A direct comparison with LSDA(+U)|| 
results clearly shows substantial discrepancy between 
theory and experiment over the whole energy range, 
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pointing to the basic inadequacy of LSDA(+U) to de- 
scribe the actual, correlated electronic structure. This 
mandates extension of LSDA to include dynamical ef- 
fects via LSDA+DMFT 

Previous theoretical work has dealt with this problem 
on two fronts. Various bandstructure calculations floL fTl| 
indicate a half-metallic state with a gap in the majority 
density-of-states (DOS). This conflicts, however, with the 
observed semiconducting resistivity up to T — 320 K, 
and so recent LDA+U work Q attempts to cure this 
malady. Moreover, the role of the B-site octahedral dis- 
tortion has been clarified in the LDA+U work, which 
suggests that the Verwey transition may be understand- 
able as an orbital order/orbital disorder transition. Ex- 
perimentally, the situation is far from clear, and theo- 
retically, one expects staggered (Neel) orbital order to 
be susceptible to strong dynamical frustration (in a way 
similar to what happens to usual antiferromagnetism on 
inverse spinels). However, the possibility of more com- 
plex orbital ordering patterns, involving more than one 
tetrahedron, is not ruled out in principle. On the second 
front, Ihle and Lorenz [l2T | explicitly considered a scenario 
involving a phase transition from a CO insulator below 
Ty to a short-range ordered, "metallic" state above Ty. 
This model gives rise to polaronic band motion above Ty 
along with polaronic hopping conductivity at higher T, in 
good agreement with observations. However, in light of 
experiments questioning the relevance (or even existence) 
of CO, a new starting point may be necessary. Motivated 
by new structural data and by the LSDA+U work || , 
an order-disorder transition in the t 2g orbital sector on 
the B-sublattice (Fe 2+ — Fe 3+ ), along with concomitant 
Jahn- Teller distortions, might be an attractive candidate 
for a consistent understanding of this system. We point 
out that study of the effects of orbital correlations/ JT 
distortions in geometric ally frustrated TM oxides is an 
enterprise in its infancy [13j ; these studies show that the 
combination of strong correlations in the d shell and ge- 
ometrical frustration can induce qualitatively new, com- 
plex behaviors. 

Summarising, the microscopic origin of carrier localisa- 
tion, and of its modification across Ty, remains an open 
issue. The clarification of this issue should go hand-in- 
hand with a microscopic understanding of the changes 
in electronic structure across the transition. Here, we 
explore this new proposal in detail. We marry the state- 
of-the-art LSDA with multi-orbital dynamical mean 
field theory [(MO)DMFT] restricted to the t 2g sector (see 
below). We choose LSDA+ (MO)DMFT because though 
LSDA captures structural effects reliably in the detailed 
one-electron dispersion(s), it generically fails to yield 
the correlated, insulating/metallic ground states charac- 
terised by appreciable dynamical SWT. 

We start with the recent LSDA density-of-states 
(DOS) [13. This describes a half- metallic ferromagnet, 
with full spin polarisation over an appreciable energy 



scale around Ep. The total DOS has dominant Fe(B) 
site contributions, but those coming from Fe{A) sites, as 
well as 0-2p states cannot be neglected either. Strictly 
speaking the Fe DOS also has contributions coming from 
d — p hybridisation. In what follows, we include only 
the Fe — t 2g DOS as the bandstructure (LDA) input 
in our calculation. The doublet E g \ orbitals lie about 
0.019 eV above the singlet A\ g orbital, giving a t 2g level 
splitting, A t2g = A = 0.019 eV. The Fe - e g states, 
split by the crystal field, lie much higher and can be 
safely neglected. Ab initio estimates yield the intra- 
orbital Hubbard U ~ 3.7 — 4.0 eV, and, with a Hund 
coupling, Jh —1.0 eV, the inter-orbital Hubbard inter- 
action, U' ~ (U - 2J H ) = 1.7 - 2.0 eV. Further, A 
acts like an external field in the orbital sector, sensitively 
controlling the occupations of the t 2g orbitals (orbital po- 
larisation) in much the same way as the magnetisation 
of a paramagnet is a function of an external magnetic 
field. In addition, the nearest neighbor Coulomb interac- 
tion, V, between n.n Fe sites is important |14| in Fe304: 
its value is not reliably known, but is estimated to be 
V ~ 0.4 eV. With this, the total many-body Hamilto- 
nian for Fea04 is, 

h = E + c/E <r<i + u ' E « 

kaber ia iab 

+ v E nf»5- J HE s ' ,s ' +A Ew- n i)' 

<i 1 j>ab iab iab 

where a, b label the t 2g orbitals. 

We use (MO)DMFT to solve this three-orbital model. 
Further, the iterated perturbation theory (IPT) is used as 
an "impurity solver". While not exact, IPT does provide 
very good agreement with QMC results for the one-band 
Hubbard model at high-T, and with exact diagonalisa- 
tion as well as very recent dynamical-DMRG results for 
the one-particle spectral functions [l^L flfij ] . The detailed 
formulation of (MO)IPT for orbitally degenerate systems 
has already been developed and used |rjj, so we do not 
repeat the equations here. Since Fea04 is half- metallic, 
we use the (MO)IPT version developed earlier by us for 
Cr02 17] to study the spectral function in the minority 
spin sector. We extend earlier work to include the ef- 
fect of V by decoupling it in the Hartree approximation, 
which is exact as d — > oo. 

For Fe304, we start with the high-T cubic struc- 
ture, with its bandstructure computed by LSDA [llj . 
In addition to generating multi-orbital Hartree shifts, 
which statically renormalise the LDA orbital energies by 
amounts depending upon U', Jh as well as by the oc- 
cupation of each t 2g orbital, the second order diagrams 
in DMFT explicitly account for the dynamical nature of 
electronic correlations, leading to drastic shifts in spec- 
tral weights and non-trivial changes in the LDA line- 
shapes. With U' = 1.7eV,V = OAeV, we derive the 
correlated spectral function in the high-T phase: the re- 
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FIG. 1: (Color online) Orbital-resolved (upper panels) and 
total (lower panel) one electron spectral functions for the two 
insulating phases of FeaCU. Notice the changes in the Ai g or- 
bital DOS; the E g \ orbitals shows small changes, demonstrat- 
ing the nearly two-fluid character of the Verwey transition. 



suit is shown in Fig. ^ The (Mott) insulating behav- 
ior is readily apparent, as is the drastic modification of 
the LSDA spectrum by multi-orbital dynamical correla- 
tions. We estimate the (Mott) gap by extrapolating the 
leading edge of the total DOS to the ui axis, obtaining 
A^ T = 0.051 eV, in good agreement with experimental 
indications. 

Motivated by the possible absence of CO in Fe304, and 
by the LSDA+U work showing up the important role of 
the JT distortion in the _B-sublattice, we adopt the fol- 
lowing strategy to study the Verwey transition. We start 
by noticing that, in a strongly correlated system, the cor- 
related values of structural parameters do not have any 
correspondence with bare LDA values anymore; in fact, 
these bare LDA parameters do not have a clear physi- 
cal meaning when correlations are strong, and are renor- 
malised in widely unanticipated ways by multi-orbital in- 
teractions. So we adopt a new strategy, where we seek 
the instability of the high-T phase as a function of the 
-B-site JT distortion. We vary the crystal field on the B 
sites in small trial steps, starting from its LDA+DMFT 
value (and not the bare LDA value) in the high-T cubic 
phase. For each trial value, we monitor the correspond- 
ing spectral functions, as well as the renormalised val- 
ues of orbital occupations and the B-site distortion from 
the correlated solution, in a way similar to that done 
for V2O3 and VO2 |l8|. Our procedure is conceptually 
very different from that used in other works where 
the two correlated phases are described by two separate 
LDA+DMFT calculations, with the corresponding LDA 
DOS for each phase as inputs into the DMFT calcula- 
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FIG. 2: (Color online) Comparison of theoretical 
LSDA+DMFT results for the total spectral function in 
the high- and low-T phases of Fes04 with experimental 
results taken from Ref.|| (for PES) and from Ref.|j (for BIS). 



tions. 

We use this strategy to investigate the Verwey tran- 
sition in magnetite. Beginning from the LSDA+DMFT 
value of A = Ea ± — Ee 1 as the initial input, we selfcon- 
sistently search for the critical value of A(= A c ) which 
destabilises the high-T phase, within LSDA+DMFT. Wc 
find, as will become clearer below, that the second so- 
lution of the DMFT equations becomes more stable for 
A c = 0.01 eV. Interestingly, (lower panel of Fig. 
the second solution also corresponds to an insulator, but 
with a larger gap O(0.1 eV): we point out that the low-T 
phase is indeed found to have a larger gap in PES 0, 0| . 
As for the high-T phase, dynamical multi-orbital corre- 
lations drastically modify the LSDA(+U) spectra in the 
low-T phase as well. 

We now compare our LSDA+DMFT results with re- 
cent published PES 0, and BIS Q measurements 
(BIS was done only for the high-T phase). Remark- 
ably, excellent quantitative agreement with the occu- 
pied as well as unoccupied parts of the spectrum is 
clearly revealed (Fig.0) over the whole energy range from 
—0.7 eV to 1.0 eV. Concomitantly, comparing the two 
LSDA+DMFT solutions, wc find that the high-T phase 
is also insulating with smaller gap (A^ T — 0.051 eV) 
than that for the low-T phase (A^ T = 0.075 eV). This 
is again in agreement with the PES spectra. 

Further light on the character of the Verwey transi- 
tion is shed by studying the evolution of the full many 
body DOS as a function of A. In Fig. |3 we show how 
the DOS depends on A within LSDA+(MO)DMFT. It 
is indeed remarkable that our results resemble the ex- 
perimental ones |5j very closely if we plausibly suppose 
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FIG. 3: (Color online) Effect of the tig level splitting, A 
on the total one electron spectral function across the Verwey 
transition in FesCU. Inset shows the effect of charge order 
(dashed line) in the low-T phase. 

that A changes with T. Up to a critical value A c of the 
crystal field splitting, all curves lie essentially on the one 
corresponding to the high-T phase. For A > A c , how- 
ever, we clearly see that all the curves seem to collapse 
onto the result for the low-T phase. Appreciable changes 
in dynamical spectral weight transfer occur across A c , in 
full quantitative accord with observations. This implies 
that the charge gap increases rapidly as A is raised above 
A c , closely tracking the resistivity jump across Ty. 

Up to this point, our calculation does not require the 
explicit inclusion of charge order (CO) within DMFT. 
This has long been thought to be essential in Fe3C>4, 
but, as argued recently (see above), even the existence 
of CO has proved to be controversial. In any case, if a 
small charge disproportionation would indeed exist, this 
would result in the addition of a minute staggered "field" 
(coming from Vnfn'j in H), or staggered chemical poten- 
tial term to the DMFT Hamiltonian. With V = 0.4 eV 
and disproportionation rj = 0.1 6], the staggered field 
term is H s = VJ(-1)X = 0.04(-l)*n? ; we expect that 
this small correction will result in miniscule changes in 
the electronic structure. However, using single site CPA 
to treat dynamical effects of CO in the low-T phase, we 
compute a small correction to the gap (at low-T) of order 
0.04(n£) ~ 0.012 (with (n\) = 0.255 for the E gl orbitals). 
From the inset of Fig. [3] the low-T gap is now estimated 
to be Ag T = 0.087 eV. The smallness of this term leads 
us to conclude that the electrostatic energy gain associ- 
ated with the CO would be too small to drive the Verwey 
transition. Moreover, the insulating behavior, as seen in 



the dc resistivity (and fitted previously by an activated 
hopping term) is also naturally accounted for: assuming 
this form, Pdc(T) ~ exp(Ag/ksT) for both phases at Ty 
with the appropriate gaps derived above, we find that the 
resistivity jump across Ty is r = e 4 ' 5 ~ 90. This is very 
close to the factor of slightly less than 100 as estimated 
experimentally |2(il |. constituting a further check on our 
results. 

In conclusion, in line with very recent experiments, we 
propose that B-site JT distortions in a correlated multi- 
orbital situation play a crucial role in understanding the 
changes in electronic structure across the Verwey tran- 
sition in Fe30<4. In the strongly correlated situation, 
small changes in the -B-site octahedral distortions with 
T give rise to large changes in dynamical spectral weight 
transfer from low- to high energies as T is raised. The 
excellent quantitative agreement with the experimental 
spectral function, as well as with the gap values in both 
phases and with the resistivity jump across Ty consti- 
tutes strong support for our new mechanism involving 
the role of multi-orbital correlations and associated struc- 
tural effects in the half-metallic situation. 
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